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In the tight-binding approximation the Harper like equation describing 
an electron in 3D crystal subject to a uniform magnetic field is obtained. 
It is supposed that the vector H can be oriented along several directions 
in the lattice. The Fermi surfaces relevant to a magnetic flux p/q = 1/2 in 
a simple cubic lattice are built. The quantization rules in magnetic fields 
slightly distinguished from p/q = 1/2 are investigated. 
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The effects of a high magnetic field on Bloch electrons are fascinating problem with 
very rich physics. This is not surprising since there are two fundamental periods in the 
problem - period of the potential and the period of the phase of the wave function - and the 
interplay of these two periods gives very interesting spectrum (Hofstadter butterfly) and 
eigenstate structure. Starting from classical works of Azbel [1], Hofstadter [2], Wannier 
[3] (see also [4-7]) this problem attracts the increasing attention. During the last decade 
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the one electron quantum states in lateral superlattices in the presence of perpendicular 
magnetic field have been studied in a series of theoretical [8-10] and experimental [11,12] 
works. But unfortunately, all works mentioned above deals only with 2D crystals subject 
to a perpendicular magnetic field. Nowadays when in VNIIEF (Sarov) the magnetic field 
up to 28 MGs is reached [13] the idea of observation of details forecasted by the theory 
in the crystals with the lattice spacing of a = 0.6 nm and more does not seem to be 
unaccessible. 

In this paper we shall consider electronic states in 3D crystals subject to ultrahigh 
magnetic field. We consider that the vector H is parallel to arbitrary translation vector of 
a lattice. In a tight-binding approximation the wave function satisfying the generalized 
Bloch-Peierls conditions is obtained and Harper like equation for different orientations of 
a magnetic field is derived. As an example, the Fermi surfaces for magnetic subbands 
relevant to rational value of a magnetic flux p/q = 1/2 are built. The quantization rules 
for magnetic flux near this rational value of flux quanta are studied. 

First of all in the tight-binding approximation let us derive Harper's equation for 
the simple cubic lattice subject to uniform magnetic field. To be specific let the basis 
vectors on the x,y plane be ai and a 2 , and the magnetic field is applied along lattice 
basis vector a 3 || Oz (|ai| = |a 2 | = |a 3 | = a). Let vector potential is be chosen in Landau 
gauge A = (0, Hx,0). Then electron wave function must satisfy the generalized Bloch 
conditions (Peierls conditions) [14] 

■0k( r ) = + qa,y + a, z + a) exp(—ik x qa) exp(—ik y a) exp(—ik z a) exp(—2iripy/a), (1) 

due to the fact that the vector potential is not a periodic function of the coordinates. Here 
p and q are mutually prime integers which define the number of magnetic flux quanta per 



two-dimensional square elementary cell 

p _ $ _ |e|if|[ai x a 2 ]| 



(2) 



q $ 2nhc 

where $o = hc/\e\ is the magnetic flux quanta. In a presence of a magnetic field the 
electron quantum states must be classified in accordance with irreducible representations 
of magnetic translation group 

&mag = nqst-i + ma 2 + Za 3 , (n, m, I are integer numbers). 

The quasimomentum k takes values in the magnetic Brillouin zone 

71 , 71 f , 1 f , 1 

< k x < — , < k v < -, <k z < -. 3 
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In the tight-binding approximation the electron wave function which satisfies the con- 
ditions (f) can be written in the form [15] 

(p (y — triaj \ 
— 2ni n\ipo{x — na,y — ma, z — la) (4) 
q a J 

where ipo(r — a n ) is the atomic function. The path-dependent geometric phase 2-7r| ( y ~™ a ^ n 
of the wave function is known to play a fundamental role in the problem. 

By substituting of (4) into Schrodinger equation and evaluating the transfer integrals 
between neighbour sights we shall obtain a system of difference equations for coefficients 
g n . The magnetic field has a non-trivial influence on the transfer integrals between neigh- 
bour lattice sights in the x direction 

A = /exp ( ± 2m^^jM* ~ «0(^(r) - U(r - a n ))^ (r - a n ) dr, (5) 



where a' n = y(n ± l)a,ma,laj, V(r) is the crystal scalar potential with periodicity a in 
three dimensions and U is the atomic potential. Calculating the transfer integrals in the 
mean-value approximation we substitute y = ma and obtain A = E where E is the 



transfer integral (5) in the absence of a magnetic field. In the y direction the transfer 
integral is 



where a„ = (na, (m ± l)a, la) and in the z direction the transfer integral is equal to E . 
As a result we obtain the well-known Harper's equation [4] 



exp(ik x a)g n+1 + exp(-i/c x a)^ n _i + 2g n cos (27m- + k y a) = e(k x , k y , k z )g n , (7) 



is the dimensionless energy measured in terms of E . The spectrum e±(k x , k y ) in Eq.(7) 
consists of q bands and depend on one single parameter p/q counting the number of 
magnetic flux quanta per unit cell (Fig.l). If flux is irrational the spectrum is a singular- 
continuum - uncountable but measure zero set of points (Cantor set). Originally the 
Eq.(7) was derived by Harper [4] using Peierls substitution [14]. 

For some values of p/q = 1/2, 1/3, 1/4 it is possible to solve the Eq.(7) analytically. 
As a result one can obtain the following expressions for dispersion laws for q magnetic 
subbands: for p/q — 1/2 (two subbands) 



B = exp 



± 2nfin) J ^ (r - a' n )(\/(r) - U(t - a n ))^ (r - a n ) dr, (6) 



where 



s{kxi ^"yi ^z^) — (^c? ^y) ^ COS Iv^Oj 



(7a) 




(8) 



for p/q — 1/3 (three subbands) 






(9) 



where (3 = \/8 — a 2 /a, a = cos3k x a + cos3k y a; and for p/q — 1/4 (four subbands) 



e { l' 2 ' 3A) (k x ,k y ) = ±\Ja ± 2^4~- sin 2 2k x a - sin 2 2k y a. 

Fig. 2 shows the Fermi surfaces in the lowest subband for p/q = 1/2 obtained using 
(7a) and (8) in the first Brillouin zone —n/qa < k x < n/qa, —n/qa <k y < Ti/qa, —n/a < 
k z < Ti j a. Three surfaces for different representative energies are plotted. For an ellipsoid- 
type surface such as (1) (|k| — > 0) the simple analytical expression for an energy spectrum 
e(k) can be easily obtained with the help of (7a) and (8). The relevant effective masses 
are determined as 

- %2 * ^ (11) 



x ' y V2E a^ z 2E a 2 ' 
The asymptotic spectrum near the point k XtV — > n/2a ((2)-type surface) is e(k x , k y , k z ) ~ 

\Jkl + k 2 — 2 cos k z . 

Now we shall discuss the magnetic quantization rules at values of a flux slightly dis- 
tinguished from p/q = 1/2: 

Hv ?,>>1 - (12) 

This corresponds to the magnetic field H = H(l/2) + AH(l/q'). As one can see from 
Fig.l near value of p/q = 1/2 the spectrum represents a system of narrow subbands 
(practically discrete levels). Here it is possible to see an equidistant spectrum, points of 
level accumulation and ranges where e± ~ \/N, (N is the level number). 

The spectrum s±(k x ,k y ) of Harper's equation near one half of flux quanta can be 
studied analytically due to the self-similarity of its structure. In the case p/q = 1/2 two 
subbands and two subsystems in the eigenvector distribution can be observed. It may 
be derived analytically from Harper's equation that every eigenvector function from a 
subsystem has an argument step An = 2. 

The magnetic flux (12) gives us the diagonal term in Harper's equation of the form 



2cos(7m + 2irn/q')g n = 2 cos 7m cos{2im/ q')g n . Now q' is a new period of the equation. 
Since cos 7m = (— l) n our system of equations (7) consists of two subsystems where only 
the sign varies. It is easy to construct an equation for each subsystem. After simple 
algebra we obtain (k x — k y — 0): 

47T 

9n+2 + 9n-2 + 2 cos — n g n = {e]_ - A)g n . (13) 

Now e]_ — 4 = e' ± is the energy in Harper's equation (7) and e± = ±\Je' ± + 4 describes 
the structure of the spectrum near one half of flux quanta and cos k z a = 0. Two signs at 
the square root correspond to two groups of energy subbands. The energy e± describes 
all parts of the spectrum near p/q = 1/2 (Fig.l). We can obtain from (13) the "square 
root" dependence at e± = 0, the clustering points e± = ±2 and Landau levels near 
e± = ±2\/2. These peculiarities of a spectrum can be obtained with the help of Onsager- 
Lifshitz quasiclassical quantization rules on the Fermi surfaces (Fig.2a,b). It is easy 
to see that magnetic quasiclassical quantization on the (l)-type surface results to the 
eqidistant energy spectrum at the bottom of the lowest magnetic subband. Accordingly, 
the magnetic Onsager-Lifshitz quantization near k z — > on (2)-type surface yields the 
"square root" spectrum. At last, it is clear that in magnetic field H || Oz there are self- 
crossing orbits on (3)-type surface. It gives the accumulation point in the spectrum e±_ 
near p/q = 1/2 (see Fig.l). It can be stressed that level spacing obtained from (13) has 
the same value that was derived using quasiclassical quantization. 

In accordance with Harper's equation (7) and (7a) the partial overlapping between q 
energy subbands is observed. When p/q = 1/2 the overlapping takes place in the energy 
interval (—2, 2). It is clear that when the electron concentrations leading to e F G (—2, 2) 
in the magnetic field H = if (1/2) + AH(l/q > ) the de Haas - van Alfen oscillations with 
different periods corresponding to the different energy level series (Fig.l) can be observed. 



Further it is necessary to make the following relevant note. The wave function (4) 
describes non-homogeneous probability distribution on the lattice sights. This is due to 
the amplitudes g n (k.) which give us non- uniform shapes. However, in real crystals due to 
Coloumb interaction the static electron density distribution must maintain its symmetry 
even in the presence of magnetic field. 

Fortunately, the spectrum of Harper's equation is g-fold degenerate since e(k x , k y + 
2nj/qa, k z ) = e(k x , k y , k z ), j = 0, 1, 2, . . . , q — 1 and a wave function with homogenious 
electron density distribution can be written as a combination of (4). Let us introduce 

<(r) = EQ( k )^^k(r), (14) 

3=0 

where ^ k (r) is from (4), Cj(k) = exp(-ik x ja), T ja ip = ip(x + ja,y, z) exp(ief /he) (/ = 
—Hyja), the operator Kj transforms ipv with k = (k x , k y , k z ) into ipv with k' = (k x , k y — 
2npj/qa, k z ). After this procedure we can write (14) in the form 

<(r) = D(k) E ex P f - 2 "~ {V ~ ^ A exp(^ka n )Vo(r-a n ), (15) 

n \ Q a J 

where D(k) = Y%=o 9n+jQ<-) does not depend on n due to the periodicity of the system. 
Now (15) describes homogeneous electron density distribution on the lattice sights and 
corresponds to the same energy as ip^. 

Let us consider now the case of a simple cubic lattice subjected to a uniform magnetic 
field oriented along the diagonal of a square in the plane (x,y). The magnetic field with 
amplitude H has now the following Cartesian coordinates: 

H = ^(1,1,0). (16) 

It is convenient to choose a new coordinate system making a rotation about the old z axis 
by 7r/4 and defining the magnetic field orientation as x 3 . This transformation is written 



as following: 

_ X 3 - Xi 
X 3 + Xi 

V = —JT' 

Z = X 2 . 

We choose a new elementary cell: it is a rectangular parallelepiped based on the new 
coordinate system vectors: a square with the side \/2a at the base in the plane (xi,x 3 ) 
and the height a. This cell is non-primitive: it consists of two atoms (additional sights are 
located at the centers of top and bottom faces). In this coordinate system the magnetic 
field is oriented along one of the elementary cell basic vectors (namely, a 3 ). According to 
the general principles in this case the classification of electron states is possible [16]. Now 
the first Brillouin zone has an area two times less with respect to the initial value. The 
magnetic field in this coordinate system is written: 

H = #(0,0,1) (17) 

and it is convenient to choose a vector potential using Landau gauge: 

A = (-Hx 2 , 0, 0). (18) 

The wavefunction in the tight-binding approximation can be written as following: 

i i \ /i \ f V x i ~~ m\p2a \ . . . 

Vk(r) = V g n {K.)< exp [2-ki — = — n) exp(2ka n )^ (r - a n ) + 

n I v Q a/V2 ' 

+ exp (2mP Xl - {m ^ 1 J 2)V ^ a n) exp(.k(a n + d))^ (r - (a n + d))). (19) 

Here ka n = k\\[2am + A; 2 an + k 3 \^2al, ipo is an isolated atomic function in the pres- 
ence of the magnetic field, the amplitudes g„(k) are the subject of further research and 
9n+q — 9n- The indexes (m, n, I) are taking all integer values independently, n numerates 
sights along x 2 axis, perpendicular to the magnetic field. The sum is taken over all lattice 



sights and two addendums in (22) correspond to the two existing crystalline subsystems. 
Each of them is a simple tetragonal lattice with the elementary cell being a rectangular 
parallelepiped with square of the side \f2~a in its base and the height a. These subsys- 
tems are displaced on the vector d with new coordinates (a/y/2, 0, a/y/2). The magnetic 
translation group is constructed by the q - time multiplication of the translation period 
along x 2 axis: x 2 — > X2 + qa. The other translations are not changed. Their periods in the 
new coordinate system are V2a which is the diagonal of a square with the side a. The 
classification of wavefunctions is possible when the number of flux quanta $/$o through 
the area a 2 /\/2 (being the minimal inter-sight area in the transversal to magnetic field 
projection) is a rational number: $/$o — p/q where p and q are mutually prime integers. 

The substitution of (19) into Schrodinger equation is performed and the transfer in- 
tegrals are calculated. We will have the equation for the energy e and the amplitudes 
#n(k): 

e^ a g n+1 + e-^n-i + 4 cos (^) cos fa^n + ^)g n = eg n . (20) 

It should be mentioned that this equation may be obtained using Peierls substitution 
k — > k — ^ with the initial spectrum at null magnetic field. In old coordinates it was 
written as e = 2(cos k x a + cos k y a + cos k z a), and in new coordinates (xi, x 2 , x 3 ) it is 

k 3 a ha 
e = 2 cos k 2 a + 4 cos cos — =. 

y/2 V2 

The equation (20) is a generalization of Harper's equation for the case of magnetic field 
parallel to a diagonal of a square which is a face of the elementary cell cube. Keeping 
a three-diagonal structure (20) has some remarkable features with respect to standart 
Harper's equation. First, (20) describes total (not "transversal" e±(k 1: k 3 )) energy as a 
function of quasimomentum k: e = e(ki, k 2 , k 3 ). Now /c 3 = const corresponds to the plane 
perpendicular to H and ki = const describes the plane parallel to magnetic field. Second, 



(20) corresponds to anisotropic Harper's equation with anisotropy ratio 

4co S (^) = A (21) 

depending on the cross-section transversal to the magnetic field which is fixed by choosing 
the value of k 3 . It is well-known that Harper's equation in the case A 7^ 2 describes 2D 
square lattice with anisotropic transfer integrals [5,6]. According to (21) the persistent 
interval for A is —4 < A < 4. In the quasiclassical limit it describes the situation where 
open-type trajectories are available. The contribution of open trajectories to the whole 
topological structure varies with the respect to a cross-section k 3 = const being chosen. 
For example, the cross-section k 3 = TT\^2/3a corresponds to A = 2 which is the case of 
isotropic Harper's equation and only close-type orbits are avaliable. This is illustrated in 
Fig. 3 where Fermi surfaces (Fig.3a,b) and spectrum (Fig.3c) are shown for p/q = 31/60. 
Then, the cross-section k 3 = ir/\/2a gives us A = which is a "full anisotropic" limit for 
Harper's equation. Here as one can see in Fig.4a,b only opened-type orbits are avaliable. 
It gives a continious energy spectrum (Fig. 4c). For remaining values of A both opened 
and closed orbits are present. In the last case the spectrum e± consists of a continious 
and quasi discrete (narrow subbands) parts. 

As a conclusion, we suggested the expression for electron wave function which satisfies 
to Peierls conditions and obtained Harper's equation from initial principles of the tight- 
binding model. The Fermi surfaces are built for different orientations of the magnetic 
field. In the following paper we plan to study the electron states in 3D crystals in the 
case of the magnetic field orientation parallel to the arbitrary translation vector of the 
lattice. The properties of de Haas - van Alfen effect in this problem will be studied. 
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Fig.l. Parts of the classical Hofstadter "butterfly" near p/q = 0, 1/3, 1/2. 

Fig. 2. The Fermi surfaces for the lower (a) and upper (b) magnetic subbands of 3D 
crystals with simple cubic lattice at p/q = 1/2, H || a 3 . (a): 1) e — —4.3; 2) e — —2.0; 3) 
e = 0.0; (b): 1) e = 4.3; 2) e = 2.0; 3) e = 0.0. 

Fig. 3. a), b) The parts of Fermi surfaces bounded by planes k 3 = ±n\^2/3a (A = 2) 
at p/q — 1/2, H || ai + a 2 for the upper subband; c) the energy spectrum £±(ki, k 2 , k 3 = 
nV2/3a) at p/q = 31/60. 

Fig. 4. a), b) The Fermi surfaces in the first Brillouin zone at p/q = 1/2 and A = 
(H || ai + a 2 ) for the upper subband; c) the energy spectrum s±(k\, k 2 , k 3 = n/aV2) at 
p/q = 31/60. 



